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Abstract 

The motion of overdamped particles in a one-dimensional spatially-periodic potential is considered. 
The potential is also randomly-fluctuating in time, due to multiplicative colored noise terms, and has a 
deterministic tilt. Numerical aimulations show two distinct parameter regimes, corresponding to free- 
running near-deterministic particles, and particles which are trapped in local minima of the potential 
with intermittent escape flights. Perturbation and asymptotic methods are developed to understand 
the drift velocity and diffusion coefficient in each parameter regime. 

1 Introduction 

The subject of this paper is the statistical characterization of the motion of particles whose (one-dimensional) 
position x{t) is governed by the Langevin equation 

dx 

= + /(^) cos{kx) + g[t) sin(fcx), (1) 

where f{t) and g(t) are random functions of time. The positive constant U is the deterministic particle 
velocity, i.e., the speed at which particles would travel in the absence of the stochastic terms, and k is 
a spatial frequency which governs how the noise sources affects the particle at position x. This equation 
may be viewed as describing overdamped one-dimensional motion in a spatially-periodic potential which 
fluctuates randomly in time, with an additional tilt. Denoting the potential as V{x^t), the equation of 
motion is 

with the potential given by 

V{x,t) = -Ux- f{t)\su\{kx) + g{t)]-cos{kx). (3) 
k ' k 

Here U is interpreted as the constant tilt of the potential. 

Equations of motion of this form arise in a number of diverse applications. In microelectronic circuit 
design applications, for example, it is important to accurately predict the effect of external noise sources 
upon the phase angle x{t) of nonlinear self-sustaining oscillators [U [2]. Under some assumptions, the 
equation of motion for x{t) has been shown [31 H] to be of the form with the tilt U representing the 
designed oscillator frequency in the absence of noise, and noise sources (voltage or current sources) being 
denoted by f{t) and g{t). The noise sources cause a diffusive drift of the phase of the oscillator, known as 
phase noise. Noise of a finite bandwidth (colored noise) may also induce a shift in the average oscillator 
frequency [4j . It has recently been demonstrated that a stochastic perturbation method may be applied 
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Figure 1: Sample trajectories from equation ([5]) shown with a single realization of the random field 
f{t)cosx + g{t)smx. The background contours and shading represent the values of the random field 
(negative field values are shaded, with dotted contours), and the lines show the solutions x{t) of equation 
([8|) with /i = 0.25 and k — 0.1 (red solid line), k = 1 (blue dot-dashed line), and k = 10 (black dashed 
line). 

to accurately predict the effects of the noise on the oscillator dynamics in the limit where U is much larger 
than the noise intensity, this being the limit of interest for engineering applications. 

Further examples of models of the form ([T]) occur in the study of weakly coupled phase oscillators [5] , 
where the noisy terms f{t) and g{t) describe fluctuations in the mean field of the oscillators [6], and in 
an agent-based model of stock market price movements [7]. For all these examples, information on the 
qualitative behavior of the first and second moments of the process x{t) described by equation ([T]) is very 
desirable. In this paper we therefore examine the mean velocity and diffusion coefficient for an ensemble 
of trajectories x{t) in all regimes of parameter space, using extensive numerical simulation, perturbation 
approaches, and appropriate non-perturbative methods. 

The random functions f{t) and g{t) appearing in equation ^ are independent, stationary, zero- mean 
Gaussian processes, which are characterized by their variance and correlation function R(t): 

{fiO)fit))^{giO)g{t))^a^R{t). (4) 

(Angle brackets are used throughout the paper to denote ensemble averages). We briefly consider the case 
where the stochastic terms f{t) and g{t) are white noise sources (so that R{t) is a delta function), but most 
of our analysis is focused on coloured noise, where the presence of a non-zero correlation timescale r for the 
noise has several nontrivial effects. Note that we do not add an independent white noise term to the right 
hand side of equation as used to model thermal diffusion in ^ for example. Neglecting such diffusion 
effects at this stage allow us to clearly demonstrate the role of the fluctuations in the potential; moreover, 
it is expected that low levels of additive white noise will not qualitatively affect the results determined 
here. 

Figure [T] provides an illustration of the effects of various parameters in the problem. These parameters 
are discussed in detail in section [3] below, but the basic features are clear from the three x{t) trajectories 
plotted in Figure[TJ If the potential tilt is sufficiently large to dominate the fluctuating part of the potential, 
then the particle moves essentially according to x{t) Ut, with small corrections due to the influence of 
the fluctuating potential. The solid red line illustrates a realization of such a case, which we class as being 
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Figure 2: Schematic of the K-fi parameter plane, showing the boundary k = max(l, /i) between the running 
regime (which hes to the left of the solid black line) and the trapping regime (to the right of the line) as 
given in equation (j23p . Symbols indicate the lines of constant-K and constant-/! along which the simulation 
results of Figures [3] and S] are collected. 

in the running regime. On the other hand, if the tilt is relatively weak compared to the fluctuating-in-time 
periodic-in-space potential, the motion of the particle is dominated by the latter. Indeed, for sufficiently 
slow fluctuations, we expect the particle to become trapped in a local minimum of the potential, with 
infrequent escapes to other nearby energy wells. The trajectory x{t) plotted as the black dashed line 
illustrates an example of this so-called trapping regime. Note that the contours and shading represent 
the gradient of the fluctuating part of the potential (shading and dotted contours denote negative values, 
while solid contours and white background indicate positive values), and the trajectory x{t) closely follows 
the position of a local minimum (zero gradient contour). We defer a full explanation until appropriate 
nondimensional parameters are introduced in section [31 but note that our main result is the theoretical 
analysis of motion in the two asymptotic parameter regimes — running and trapping — introduced here, and 
depicted schematically on the parameter plane of Figure [H 

The remainder of this paper is structured as follows. In section [2] we examine the white-noise limit 
for the random processes f{t) and g{t) in equation ([3]). More interesting results emerge when the noise 
processes have a non-zero correlation time r; the analysis of such colored noise cases begins with the 
introduction of nondimensional parameters in section [3] Methods of numerical simulation are detailed in 
section |4] and theoretical analysis of the two main parameter regimes is undertaken in sections [5] and [G] 
Section [7] concludes with a discussion of results. 

2 White noise 

If the noise sources /(t) and g{t) in equation ([T]) are white (i.e., R{t) — 6{t) in equation (|4])), then the 
probability density function P{x,t) giving the particle concentration at position x at time t is the solution 
of a Fokker-Planck equation [5]. Some care is needed when interpreting this white- noise limit, as the 
noise sources appear in equation ([T]) in a multiplicative fashion. It is typical in such problems that the 
Ito and Stratonovich interpretations of the white noise terms in equation ([3]) lead to different results. 
Rather remarkably however, in this particular example the same Fokker-Planck equation is found using 
both interpretations, as the symmetry of the two noise sources causes the so-called spurious drift term to 
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vanish. Thus either the Ito or Stratonovich interpretation of the white noise hmit of equation ((T|) leads to 
the following equation for P{x,t): 

It is straightforward to determine from this equation that the mean value of x{t) grows at a constant rate 
U; indeed, multiplying equation ([5]) by x and integrating over space yields the exact relation 

(xit)) = {x{0)) + Ut, (6) 

so that the mean velocity is just U. Similarly, the diffusion coefficient of x equals the coefficient of the final 
term of equation ([5]), i.e. Note this coefficient shows no dependence on the value of the potential 

tilt parameter U. As we show below, qualitatively very different results (including diffusion coefficients 
which depend on U) are found when the noise sources have a non-zero correlation timescale. We therefore 
proceed to consider the case of colored noise in the next section. 



3 Colored noise and nondimensional parameters 

A physically more realistic model of random fluctuations than the inflnitely-fast variation of white noise 
is provided by colored noise, for which the correlation function R{t) defined in equation ^ decays from 
its peak value R{0) = 1 to zero over some characteristic timescale r. One example of such a correlation 
function is 

R{t) = e-^. (7) 

This correlation function is particularly convenient for numerical simulation, and also has the property 
R'{0) = which is important for some of our theoretical analysis below. This smoothness of R at zero 
argument is related to the finiteness of the variance of the processes f'{t) and g'{t), but it is not always 
guaranteed: if f{t) and g{t) are Ornstein-Uhlenbeck processes, for example, then this smoothness condition 
is violated. Nevertheless, for the purposes of this initial examination of the properties of the equation ([2]) 
we assume the smoothness property when needed, and use the correlation function ([7]) in all numerical 
simulations. 

Having determined a characteristic timescale r, it proves convenient to introduce dimensionless variables 
(temporarily denoted by tildes) by measuring time in units of t and space in units of 1/k. The variance 
of the noise sources has dimensions of (length/time)^, while the tilt U has dimensions length/time. The 
dimensionless equation of motion for x{t) may therefore be written as 



dx 



/(f) cos X + g{i) sin x , (8) 



where the dimensionless tilt is 
and the Kubo number is defined as 



At = Ukr (9) 

k — akr. (10) 

Note the dimensionless noise sources / and g in equation ^ have unit variance and unit (dimensionless) 
correlation time. The physical meaning of the Kubo number can be understood best in the case of zero 
tilt, i.e., with U ~ Q in equation ([T]). Then particles moving in the fiuctuating potential will travel a 
(dimensional) distance of order ar during the correlation time of the noise processes; the Kubo number is 
the ratio of this lengthscale to the reciprocal spatial frequency k~^ of the potential. Thus the h ^ limit, 
for example, corresponds to a potential which fiuctuates rapidly in time, so that a particle moving in it 
travels a very short distance (compared to the period of the potential) during the time in which values of 



4 




Figure 3: Normalized mean velocity (a) and diffusion coefficient D (b) as functions of n/fi for fixed 
values of k. In panel (a), solid lines show equation (|32p for the mean velocity in the running regime, while 
the dot-dashed line is the trapping regime model result ((54)) . In panel (b), the solid line is the leading-order 
perturbation result (|38|) for D in the running regime, while the dashed lines include the order-K^ effects 
given in equation (|44p . The dot-dashed line is the trapping regime value of D, given by equation (|48p . 



f{t) and g{t) change significantly. On the other hand, taking k — > oo corresponds to a frozen (quenched) 
limit, where particles feel the spatial variation of the potential much sooner than the temporal variation. 

For the remainder of the paper we concentrate on the dimensionless form of the equation ([8]) and 
henceforth drop the tilde notation. We focus particularly on the range of behaviors given by varying the 
two positive parameters fj, and k. Results will be presented for the large-time mean velocity v and the 
diffusion constant D of an ensemble of particles. The mean velocity is defined as 

V = hm ^ {x{t)) , (11) 

or equivalently by the large-time asymptotic form of the mean position: 

{x{t)) ^vt as t~*(x,. (12) 

Assuming diffusive behavior over sufficiently long timescales, the diffusion constant is defined by the large- 
time asymptotic relation 

var(a;) ~ 2Dt as t ^ oo, (13) 
where var(a;) is the variance of the distribution of particle positions at time t: 

V8ir{x) = (x^t)) - (x{t)f . (14) 

Alternatively, we may use the limit definition 

D = — lim — var(x). (15) 
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Figure 4: Normalized mean velocity v/fi (a) and diffusion coefScient D (b) as functions oi n/fj, for fixed 
values of /i. In panel (a), solid lines show equation (|32p for the mean velocity in the running regime, 
while the dot-dashed lines are the trapping regime model result ([54|) . In panel (b), the solid lines are the 
leading-order perturbation results (I38p for D in the running regime, while the dashed lines include the 



order-K effects given in equation (|44|) . The dot-dashed line is the trapping regime value of D, given by 
equation 



4 Numerical simulations 

Gaussian random functions may be constructed using a combination of a large number N of Fourier modes, 
as follows: 

1 ^ 

f{t) = ^ a„ cos(w„t) -f 6„ sin(w„i). (16) 

n=l 

The amplitudes a„ and 6„ are random numbers from a Gaussian distribution of mean zero and unit 
variance. The random frequencies aj„ are from a distribution shaped as the Fourier transform of Kif) 
— for the correlation function ([7]) this means the tj„ are also Gaussian distributed, with mean zero and 
(nondimensional) variance one. The function ]{t) constructed in this way is Gaussian in the limit N oo 
[TT| I12j: we use N = 100 in the numerical simulations reported here. 

In each realization, the independent noise functions f{t) and g{t) are generated as above. The ordinary 
differential equation ([5]) is then solved numerically, with initial condition x{0) = 0, to determine x(t). By 
averaging over a large ensemble (typically 10^) of realizations, we determine the time-dependent moments 
{x{t)) and (a;^(t)), and the variance var(a;). A linear fit to the average position and variance at large t (we 
use nondimensional end-times of t = 10 to t = 50 for this fitting) gives, respectively, the mean velocity v 
and the diffusion constant D from the slope of the fitted lines, in accord with the definitions (fT2|) and (I13p . 

Results from these numerical simulation are plotted as symbols in Figures [3] and ID Note the same 
symbol types are used as in the schematic parameter plane of Figure [2l with the constant-K results shown 
in Figure [SI and the constant-/! results in Figure ID We choose to plot all results as functions of k/^ in 
order to highlight several important scalings. Panel (a) of each figure shows the normalized mean velocity 
v/fi, while the diffusion coefficient D appears in the panel (b). Note all scales are logarithmic. It appears 
that two distinct regimes of simple scalings arise: broadly speaking these regimes appear to correspond to 
k/ii being respectively much smaller than, and much larger than, unity. However, as we show in sections 
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[5] and [6] below the detailed picture is somewhat more complicated, with the border between the running 
regime (section [S]) and the trapping regime (section [5]) having the shape shown in the schematic Figure 
[51 Analytical results are derived below for limiting behaviour in each regime and these are shown as the 
various curves in Figures [3] and [31 

5 Running regime and perturbation method 

When the mean velocity /i is sufficiently large, or the noise parameter n is sufficiently small, particles 
moving in the fluctuating potential according to equation ([8]) may be treated as slightly perturbed from 
their deterministic paths. In this parameter regime (to be defined in detail below) particles typically move 
with a speed close to /i and move swiftly over the potential landscape (i.e., covering many spatial periods 
during the fluctuation timescale of the potential); we therefore call this the running regime of parameters. 
This regime is of particular interest to designers of oscillator circuits and some of the results of this section 
have been previously derived (using different methods) in [J . 

To study this regime more carefully, it is helpful to move with a reference frame traveling at the 
deterministic mean velocity /x. By introducing the new variable y{t) defined as 

y{t) = x{t) - nt, (17) 

equation ([8]) may be rewritten as 

y = n [fit) cos (y + fit) + g{t) sin [y + fit)] . (18) 

We combine the noise sources and time-oscillatory terms to create two new stationary Gaussian random 
processes F{t) and G{t), defined by 

F{i) — f (t) cos ijt + g{t) sin fit, 

G{t) = ~ f (t) sin fit + g{t) cos fit. (19) 

From these definitions, it is straightforward to verify that F{t) and G{t) each have mean zero, and corre- 
lation function 

(F{t)F{0)) = (G'(t)G(O)) = R{t) cos fit. (20) 

However, unlike f{t) and g(t), the new source terms F{t) and G{t) are not independent, since their cross- 
correlation is given by 

(F(t)G(O)) = - (G(t)F(O)) = R{t) sin fit. (21) 
In terms of the new noise processes, the equation of motion (1181) is 

y = K[F{t)cosy + G{t)siny]. (22) 

Next, we assume a regular perturbation expansion for y(t) in powers of the parameter k. The results 
of such an approach will be useful if the distance traveled by the particle in the characteristic time for 
the potential to decorrelate is significantly less than the spatial period of the potential. Provided this 
condition holds, the errors in approximating cosy and siny by their series expansions about y = do not 
accumulate over time. In terms of the dimensionless parameters introduced in section [31 the spatial period 
of the potential is of order one. The decorrelation time of the potential felt by the particle in the moving 
reference frame depends on fi; from (|20|) we can estimate this as min(l, 1/fi). The distance traveled by the 
particle during this time is of order Kmin(l, 1/fi) and the condition for the perturbation expansion to hold 
is therefore expressible as Kmin(l, 1/fi) <ti 1, or equivalently as 

K<max(l,/i). (23) 
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We adopt this condition as defining the running regime of parameter values. 

We proceed to assume a regular perturbation expansion for y(t) for the k limit of the form 

y{t) = Kyi{t) + K^y2{t) + K^ysit) + ..., (24) 

where each yi(t) is a stochastic process. Substituting this into equation (|22p . expanding the trigonometric 
terms, and gathering powers of k yields the following system of equations for the unknown functions yi{t), 
y2{t), and ysit): 

m = F{t) (25) 
h = G{t)yi (26) 

2)3 = -^yl + G{t)y,. (27) 

5.1 Mean velocity in running regime 

The solution of (j25p with initial condition j/(0) = is 

yi{t) - / F(h)dh (28) 







t 

G{t)F{ti)dti. (29) 



and substitution into ([26|) yields 

dy2{t) 
dt 

Ensemble-averaging equations (|28l) and (l29|) leads to 

(yiW> = 0, 

d f* 
dt Jo 

= - [ R{s)sm{fis)ds, (30) 

where the substitution s ^ t — ti has been used. This result implies that in the long-time limit t ^ oo, we 
obtain a finite mean value for d {y) /dt at leading order ac^, given by 

d P°° 

lim — (y) = -K^ / R{s) sin(^s) ds + 0{k^). (31) 

t^OO dt Jq 

The mean velocity defined in equation (jlip is therefore given to leading order (as k — > 0) by 

d 

v= lim — (x)=/i-K^ / R{s)fiin{ij.s)ds + 0{K'^). (32) 
t^oo dt Jq 

The normalized mean velocity v //z given by this result is plotted as a solid line in panels (a) of Figures [3] 
and [4] It matches numerical results well within the running regime defined by equation (|23p . 

5.1.1 Large- /i asymptotics of mean velocity 

To investigate the asymptotic behavior of |32|l at large fj, we integrate by parts: 

/ R{s) sm{^is) ds = - R{t)-cos{fit) +-/ R' (s) cos{fis) ds. (33) 
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For sufficiently smooth correlation functions R{t), this integration by parts may be repeated in an iterative 
fashion, to obtain an asymptotic series is inverse powers of /i 13J. To leading order in we obtain 

/ R(s)sm(fis)ds^— + o(-^] as ^ ^ c» (34) 
Jo M VM / 

and we conclude that deviations of the mean velocity from ^ are of order in this regime of parameter 
values. 

5.2 Diffusion coefficient in running regime 

The variance of y{t) (which equals the variance of x{t)) is given to leading order by 

ve.Ty = K^yl{t))+0{K^), (35) 

and from (pS)) we obtain 

{yl{t))= [ dh f dt2R{tl-t2)cOS^l{h-t2). (36) 

Jo Jo 

The slope of var y(t) is given by its derivative with respect to t: 

d /■* 

—var y = 2k / R{s) cosies) ds + 0{k*) (37) 
dt Jq 

and as t — > oo this reaches the finite limit giving D from equation psp : 

/>oo 

D = K^l R{s)cos{fis)ds + 0{K'^). (38) 



This leading-order result for D is plotted with solid lines in panel (b) of Figures [3] and ID While it matches 
numerical results well for low k and low fi cases, it is noticeable that even well within the running regime 
(with K <C /i) it is very inaccurate (not even appearing on the scale of the Figures) for k = 10 (Fig. Eljb)) 
and for ^ = 10 (Fig. IHJb)). The reason for this loss of accuracy is examined in detail in the following 
section. 

5.2.1 Large- /i asymptotics of diffusion coefficient 

We consider the asymptotic behavior of the diffusion constant at large values of fi. For correlation functions 
R{t) which are sufficiently smooth at t = 0, the order integral in equation ([55]) can decrease rapidly 
with increasing /i. For example, using the correlation function of equation ([7]) in equation (j38p leads to the 
prediction that the diffusion constant is exponentially small at large /i values: 



D = -e^— +0{k^) a.s n ^ oo. (39) 

Because the first non- vanishing term in the perturbation expansion for D may be exponentially small, we 
proceed to calculate to higher orders. The contribution to var y at order is given by 

(vl) + 2 {ym) - {y2f , (40) 
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and after extensive algebraic manipulation this can be shown to yield a contribution to the diffusion 
constant D of 



D2 = K dsi / ds2 / ds 



oc 



{i?(si + S2 + 53)^(52) [—3 cos/i(si + S2 + S3) COS /iS2 + sin /i(si + S2 + S3) sin ^82] + (41) 
R{si + S2)R{s2 + S3) [-2cos/^(si + S2)cos^(s2 + S3) + 2sin^(si + S2)sin^(s2 + S3)]} . (42) 

The large-// asymptotic behavior of this term may be calculated (in a similar fashion to the mean velocity 
case above) as 

^4 roc 

D2-— R\t)dt as 00. (43) 

Jo 

For the correlation function of equation ^ we therefore have the following small-K, large-/i asymptotic 
behavior: 



^ " '^V^'^"'^ ^ /^"^ Ai->oo, K^O. (44) 

This corrected form for D is plotted with dashed lines in Fig. ^h) and Fig. Iljb) and shows excellent 
agreement with numerics within the running regime. We conclude that the deviation from the leading- 
order result (|38p noted above is due to the order-K^ effects dominating the order-K^ term when /i is 
sufficiently large. 



5.3 Summary of running regime 

In brief, condition (j23p defines the running regime of parameter values and within this regime the mean 
velocity and diffusion coefficient may be determined using a stochastic perturbation expansion. To order 
the mean velocity v is given in equation ((32)l and the diffusion coefficient in equation p8|) . Asymptotic 
analysis for large ^ reveals that the velocity scales as /i — as /i ^ cx), while D can have an order 
term which is exponentially small, see equation (j39p . In this case the order-K^ contribution of equation 
(021) will dominate. 



6 Trapping regime 

In section [S]wc analyzed particle motion under the condition that k <^ max(l,/i), which corresponds to 
the running regime. At the other extreme, i.e. when n ^ max(l,/i), the dynamics of particles may be 
understood in terms of trapping at low-velocity positions, with intermittent releases and flights between 
trapping points. To analyze the dynamics in this trapping regime, it is useful to rewrite equation ([8]) as 

dx 

— = fi + Kr(t) sin (x + (t>(t)) . (45) 
dt 

The phase (j)(t) and the amplitude r{t) are each random functions of time, and are given in terms of the 
original noise processes f{t) and g{t) as 

cf>{t) = tan-i r{t) - ^JW+gW; (46) 

conversely, the noise sources may be expressed as f{t) = rsin(j) and g = rcoscj)- 
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6.1 Diffusion coefficient in trapping regime 

Taking the limits — > and k — > oo in equation (j45p allows us to derive some useful results. In this 
limit, the particle trajectories x{t) are locked to the phase variable: x{t) — —4>{t) (up to the addition of 
multiples of 27r), and in particular the velocity of the particle is given explicitly as a function of time by 
x{t) = —4>{t). The phase velocity 0(i) in systems of this type has been extensively studied [TH [6], and 
known results may be applied to find the diffusion coefficient in this limit [7j. The mean value of </> is zero, 
and the phase velocity correlation fmiction may be expressed in terms of the correlation function R.{t) of 
the noise sources as 

m = (mkt)) = [R{t)R"{t) R'{tf] log [1 - R{tf] . (47) 

In this trapping limit the effective diffusion coefficient for the particles may therefore be approximated by 
setting x{t) = — (?!>(i) and proceeding as follows: 



D = llun^{x^t)) = lunUimt)) 

2 t^oo at t^co \ I 



lim ( / (j){ti)(t){t) dt 







L{ti)dti (48) 







Given the noise correlation function, this integral may be calculated numerically after applying equation 
(|T7)) : for the case R{t) = exp(— this yields the value D = 1.158. This value is plotted with a dot- 
dashed line in Figures Elb) and (Hb) and agrees with the numerical results for parameters well within the 
trapping regime k ^ max(l, /i). 

6.2 Mean velocity in trapping regime 

In the strong trapping limit of /i — > the mean velocity of the particles is necessarily zero. Numerical 
results (e.g. Fig. HJ^a)) indicate that if fi is non-zero but small then the mean velocity scales linearly with 
(i/k, and we are motivated to study the behavior of equation (|45p as a small-/i perturbation about the 
strong-trapping case considered above. As we shall see, it is a challenging problem to describe the dynamics 
in this case, because of short (in time) but large (in space) intermittent escapes of the particle from the 
strong-trapping solution. As a consequence, we will revert to a simpler toy problem in an attempt to 
understand the observed numerical scaling of the mean velocity in this regime. We begin by giving the 
exact equation for the deviation z{t) — x(t) + (j)(t) of the particle position x{t) from its strong-trapping 
limit ~(f>(t): 

dz 

— ^ fi + (f){t) + Kr{t) sin z. (49) 

Analysis of this equation is significantly complicated by the dynamics of the non-Gaussian phase velocity 
(j){t). For example, noting that the phase velocity may be expressed as 

(p — — tan — 
dt g 

gf - fg 

(50) 



f cos (f> — g sin ( 
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and that f{t) and g{t) are independent Gaussian processes, we see that the one-time distribution of (f)(t) is 
the same as that of a(t)/r{t), where a{t) is a zero-mean and unit variance Gaussian random function and 
r{t) is the ampHtude given by equation (^5]) above. This fact may be used to find the distribution of (f>{t) 
in closed form [7j, but here we use it to develop an analytically-tractable toy model which demonstrates 
similar behavior to the full dynamics. 

We consider splitting time into intervals of length T, and in each interval replacing the time-varying 
functions r{t) and (f>(t) with random constants chosen from the appropriate one-time distributions. This 
frozen-noise approximation thus neglects all variation of the phase and amplitude variables on timescales 
shorter than T, but it allows the dynamics of the corresponding z{t) to be found be solving the constant- 
coefficient equation 

dz (7 , 

— = n -\ hKrsinz (51) 

at r 



within each T-interval. Here r is given by r = y\P+g^ and is approximated by cr/r, with /, g, and 
a being independent Gaussian random numbers with mean zero and unit variance, chosen anew for each 
new time interval. 

We focus on the mean velocity which is equal to {z{t)) since </> has zero mean. The ensemble- 

averaged z-velocity within the toy model is calculated by solving (|5ip within each T-interval, and then 
averaging over the possible values of r and a. There are essentially two cases to consider. If the random 
numbers a and r satisfy the inequality 

fi+-Y~K^r^<0 (52) 
r / 

then a steady state solution of ([5T|) exists for z(t), and so after a transient the particle reaches this steady 
state and the velocity decays to zero. For sufficiently long time intervals T, we can neglect the transient and 
so consider the effective velocity to be zero in cases satisfying ([52| . If the inequality is not satisfied, however, 
then no steady state solution of equation (|5ip exists, and these running solutions have a time-averaged 
velocity calculated (see the appendix) to be 



r 



To determine the ensemble mean velocity it remains only to average these running velocities over the 
distribution of r and a values satisfying the inequality (|52p . After some manipulation (see appendix) it is 
possible to show that the resulting mean velocity is linear in /i, taking the form 

(.).iy|^ + 0(.= ) .s^^O^ (54) 

Using this form for the mean velocity gives (after normalizing by n) the dot-dashed curves in Figures [3)Ja) 
and|4ja). It is apparent from Fig. |4l^a) that the linear dependence on /u/k in (|54p qualitatively matches the 
numerical results, but the prefactor is too small by a factor of approximately 2. Given the approximations 
made in the simplification to the toy problem used here, this level of agreement is considered rather 
satisfactory. 



7 Discussion 

To summarize: we have examined the behavior of first and second moments of the solution x{t) of the 
equation ([1]) (nondimensionalized to equation ([S])) with independent random forcing terms f{t) and g(t). 
Numerical simulation results are supported by asymptotic analysis in the two main regimes of the K-fi 
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parameter plane, termed the rumiing regime (section [S]) and the trapping regime (section [6]) — see also 
Figure [Hand definition 

The smoothness of the noise correlation function R{t) at t = was assumed at the beginning of section 
[31 If the smoothness property R'{0) = is violated by, for example, using Ornstein-Uhlenbeck processes 
for f{t) and g{t), then the analysis of the trapping regime which relies on f{t) and g{t) being finite- variance 
processes is not valid. The perturbation analysis in the running regime is unchanged, but note the leading- 
order diffusion term typically scales as k^/h^ for small k and large /i if i? is non-smooth and so the 
order-K^ term P5)) will not dominate as it does when R is smooth. 

It is interesting to compare the results for our model with related work on noisy transport in periodic 
potentials. The non-monotone dependence of the diffusion coefficient on k/^, seen in Fig.[3]Jb) for k = 10 
and in Fig. [JJ^b) for fi = 10, is reminiscent of Figure 1 of [8]. The authors of [8] consider overdamped 
motion in a tilted, spatially-periodic potential which (unlike our randomly-fluctuating case) is static in 
time, but to which white noise effects (not included by us) are added. They find the diffusion coefficient 
in their model is greatly increased at a critical value of the tilt parameter corresponding to the onset of 
deterministically running solutions. We find the peak value of D in our model occurs for k/ii ratios of 
order unity (actually approximately 0.5 to 0.7), which corresponds to the border between the running and 
trapping regimes at high values of n or fi. Because particles may be either running or trapped in any given 
realization, the overall variance of x{t) across the ensemble is particularly large for n/fj, fa 1, leading to a 
diffusion coefficient which exceeds that found in cases where each single realization is more typical of the 
ensemble average. 

Our perturbation expansion method for the running regime bears some resemblance to that used in 
[l5] to examine Stokes drift in a fluid due to the presence of periodic waves. Substitution of the Fourier 
expansions (fT6|) for f{t) and g{t) in equation ([8]) allows the equation of motion to be written as 

^ = A* + -y= ^ An C0S{X - UJnt) + Bn Sm{x - LUnt), (55) 
* n—1 

with the sum admitting the interpretation as a superposition of waves of various amplitudes (and directions) 
and random frequencies w„. Our result (|32p may therefore be compared (by taking = 0) to the drift 
velocity determined by a perturbation exapansion in [15j . Unlike [15] . however, the randomness of our 
waves means there is no preferred wave direction, and so we find a zero drift velocity (we also do not 
consider the additive white noise diffusion effects used in [15] and [l6j). The directional symmetry of the 
random wave field is broken when is non-zero and in this case the deviation of v from the deterministic 
velocity fi is of order (from equation (|32)l ). which (like in [l5j) is the square of the wave amplitude. It 
would be of considerable interest to extend our methods to include the effects on the particle trajectories 
of additive white noise terms as in [8j, [15] . [16] . in addition to the multiplicative colored noise terms 
considered here. 

Finally, we note the possibility of using spatially-periodic potentials for sorting of microparticles, which 
has attracted considerable recent attention [T7] , [TH] , [IH] , [10] ■ If experimental periodic potentials can also 
be made to fluctuate randomly in time (perhaps by modulation of the laser intensities generating optical 
lattices as in [17j . |18| ) then the results of this paper would provide a useful basis for predicting mean 
velocities and spreads of particles being sorted within the device. 
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Appendix: Mean velocity in toy model for trapping regime 



In section 16.21 we introduced a simplified toy model for the trapping regime within which we can gain 
some analytical insight into the behavior of the mean velocity. Recall equation (|5ip generates a non-zero 
contribution to the mean velocity (i) in a given T-interval if the parameters r = \/ + and a do not 
satisfy inequality (j52p — here /, g and a are unit variance, mean zero normal random numbers in each 
T-interval. Indeed, in this case equation (|5ip immediately yields an integral solution for the time taken to 
travel the spatial period 2it: 

— ^- — = / dt, (56) 

/i + — + Kr sm z Jq 

which gives the equation (|53p in the main text for the running velocity Vmn in that particular T-interval 
(assuming T is sufficiently large). In all cases the sign of Vmn is the same as the sign of the quantity /U -h ^. 

Next, we consider the averaging of (|53p over the possible values of r and a. The values of a are chosen 
from the range (— cxd,cxd) with density e~°^ /^/•\/27r, while the values of r are from the range [0,cxd) with 
density re~^ 1"^ . When averaging over all values of r and a satisfying the inequality 

Ai+^)^-KV2>0 (57) 
we must consider three regions of the r-tr parameter plane, as follows: 



Region 1: cr = to oo, with r = to (/i + -^//i^ + 4k(t)/2k. In this region Wjun is positive. 



• Region 2: a ~ —[i^j^n to 0, with r ^/y?+Ana)/2n to (/i -f \/ \J? ^ Ana) /{2k). In this region 
Viun is positive. 

• Region 3: a — — oo to 0, with r = to (— /i + fi'^ — Ana) /2k. In this region Vi^n is negative. 

Because we are interested only in the leading-order behavior of (i) at small deterministic velocity /i, we 
focus on the linear term in an expansion in small fi: 



+ 0{fi'). (58) 



It is easy to show that Region 2 does not contribute to the mean velocity at this order, while Regions 1 
and 3 each contribute an amount equal to 



27r Jo Jo d^i 



dr da 

M=0 



V2tt Jo Jq \fa 
After the change of variable w = kt^ / a we obtain 



' ' e 2-rp~~ - drda (59) 



2k Jq Jq a/T 

^ [tt f 1 



ere 2 e « : dw da 

, - , , as K ^ oo, (60) 
4k V 2 V « 

where the final expression has been obtained by taking k ^ oo in the integrand. Adding the equal 
contributions from Regions 1 and 3 then gives the asymptotic expression of equation (j54p of the text for 
the averaged velocity. 
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